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Abstract. The string equation of type (2, 2g + 1) may be thought of as a higher order 
analogue of the first Painlevc equation that corresponds to the case of g = 1. For g > 1, this 
equation is accompanied with a finite set of commuting isomonodromic deformations, and 
they altogether form a hierarchy called the PI hierarchy This hierarchy gives an isomonod- 
romic analogue of the well known Mumford system. The Hamiltonian structure of the Lax 
equations can be formulated by the same Poisson structure as the Mumford system. A set 
of Darboux coordinates, which have been used for the Mumford system, can be introduced 
in this hierarchy as well. The equations of motion in these Darboux coordinates turn out to 
take a Hamiltonian form, but the Hamiltonians are different from the Hamiltonians of the 
Lax equations (except for the lowest one that corresponds to the string equation itself). 

Key words: Painlevc equations; KdV hierarchy; isomonodromic deformations; Hamiltonian 
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1 Introduction 

The so called 'string equations' were introduced in the discovery of an exact solution of two- 
dimensional quantum gravity [H EJ [3]. Since it was obvious that these equations are closely 
related to equations of the Painleve and KdV type, this breakthrough in string theory soon 
yielded a number of studies from the point of view of integrable systems [H El [6j [8] . 

The string equations are classified by a pair (q,p) of coprime positive integers. The simplest 
case of (q,p) = (2, 3) is nothing but the first Painleve equation 

1 3 2 

2 u xx + + x = 0, 

and the equations of type (2,p) for p = 5, 7, . . . may be thought of as higher order analogues 
thereof. Unlike the case of type (2, 3), these higher order PI equations are accompanied with a 
finite number of commuting flows, which altogether form a kind of finite-dimensional 'hierarchy'. 
This hierarchy is referred to as 'the PI hierarchy' in this paper. ('PI' stands for the first Painleve 
equation) . 

The PI hierarchy can be characterized as a reduction of the KdV (or KP) hierarchy. This is 
also the case for the string equations of all types. The role of the string equation in this reduction 
resembles that of the equation of commuting pair of differential operators [ID]. It is well known 
that the equation of commuting pairs, also called 'the stationary Lax equation', characterizes 
algebro- geometric solutions of the KP hierarchy I12| . The reduction by the string equation, 
however, is drastically different in its nature. Namely, whereas the commuting pair equation 
imposes translational symmetries to the KP hierarchy, the string equation is related to Virasoro 
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(and even larger Wi+oo) symmetries [6]. In the Lax formalism of the KP hierarchy |13^ 114} [T5]. 
the latter symmetries are realized by the Orlov-Schulman operator [TB], which turns out to be 
an extremely useful tool for formulating the string equation and the accompanied commuting 
flows [9]. 

The string equation and the accompanied commuting flows may be viewed as a system of 
isomonodromic deformations. This is achieved by reformulating the equations as Lax equations 
of a polynomial L-matrix (a 2 x 2 matrix in the case of the PI hierarchy). This Lax formalism 
may be compared with the Lax formalism of the well known Mumford system |17j . Both systems 
have substantially the same 2x2 matrix L-matrix, which is denoted by V(X) in this paper. A is 
a spectral parameter on which V(X) depends polynomially. The difference of these systems lies 
in the structure of the Lax equations. The Lax equations of the Mumford system take such 
a form as 

d t V(X) = [U(X),V(X)], 

where U(X) is also a 2 x 2 matrix of polynomials in A. (Note that we show just one of the Lax 
equations representatively.) Obviously, this is an isospectral system. On the other hand, the 
corresponding Lax equations of the PI hierarchy have an extra term on the right hand side: 

d t v{x) = [c/(A), v(x)} + u'(x), u\x) = d x u(x). 

The extra term U'(X) breaks isospectrality. This is actually a common feature of Lax equations 
that describe isomonodromic deformations. 

We are concerned with the Hamiltonian structure of this kind of isomonodromic systems. As 
it turns in this paper, the PI hierarchy exhibits some new aspects of this issue. Let us briefly 
show an outline. 

The Hamiltonian structure of the Mumford system is more or less well known [18^ [19] . The 
Poisson brackets of the matrix elements of the L-matrix take the form of 'generalized linear brac- 
kets' [20j. (Actually, this system has a multi-Hamiltonian structure [21\ I22j. but this is beyond 
the scope of this paper.) The Lax equations can be thereby expressed in the Hamiltonian form 

d t V(X) = {V(X),H}. 

Since the Lax equations of the PI hierarchy have substantially the same L-matrix as the 
Mumford system, we can borrow its Poisson structure. In fact, the role of the Poisson structure 
is simply to give an identity of the form 

[U(X),V(X)] = {V(X),H}. 

We can thus rewrite the Lax equations as 

d t V(X) = {V(X),H} + U'(X), 

leaving the extra term intact. This is a usual understanding of the Hamiltonian structure of 
isomonodromic systems such as the Schlesinger system (see, e.g., [231 Appendix 5]). This naive 
prescription, however, leads to a difficulty when we consider a set of Darboux coordinates called 
'spectral Darboux coordinates' and attempt to rewrite the Lax equations to a Hamiltonian 
system in these coordinates. 

The notion of spectral Darboux coordinates originates in the pioneering work of Flaschka and 
McLaughlin [23], which was later reformulated by Novikov and Veselov [25] in a more general 
form. As regards the Mumford system, this notion lies in the heart of the classical algebro- 
geometric approach [T7]. The Montreal group [26] applied these coordinates to separation of 
variables of various isospectral systems with a rational L-matrix. Their idea was generalized 
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by Sklyanin [27] to a wide range of integrable systems including quantum integrable systems. 
On the other hand, spectral Darboux coordinates were also applied to isomonodromic systems 

ESI EDI EH- 

We can consider spectral Darboux coordinates for the PI hierarchy in exactly the same way 
as the case of the Mumford system. It will be then natural to attempt to derive equations of 
motions in those Darboux coordinates. Naive expectation will be that those equations of motion 
become a Hamiltonian system with the same Hamiltonian H as the Lax equation. This, however, 
turns out to be wrong (except for the lowest part of the hierarchy, namely, the string equation 
itself). The fact is that the extra term U'(X) in the Lax equation gives rise to extra terms in 
the equations of motions in the Darboux coordinates. Thus, not only the naive expectation is 
negated, it is also not evident whether those equations of motion take a Hamiltonian form with 
a suitable Hamiltonian. This is the aforementioned difficulty. 

The goal of this paper is to show that those equations of motion are indeed a Hamiltonian 
system. As it turns out, the correct Hamiltonian K can be obtained by adding a correction term 
AH to H as 



This is a main conclusion of our results. It is interesting that the Hamiltonian of the lowest flow 
of the hierarchy (namely, the string equation itself) is free from the correction term. 

Let us mention that this kind of correction terms take place in some other isomonodromic sys- 
tems as well. An example is the Gamier system (so named and) studied by Okamoto [32J. The 
Gamier system is a multi-dimensional generalization of the Painleve equations (in particular, the 
sixth Painleve equation), and has two different interpretations as isomonodromic deformations. 
One is based on a second order Fuchsian equation. Another interpretation is the 2x2 Schlesinger 
system, from which the Gamier system can be derived as a Hamiltonian system for a special 
set of Darboux coordinates. It is easy to see that these Darboux coordinates are nothing but 
spectral Darboux coordinates in the aforementioned sense \28\ I29j. and that the Hamiltonians 
are the Hamiltonians of the Schlesinger system [231 Appendix 5] plus correction terms. These 
observations on the Gamier system have been generalized by Dubrovin and Mazzocco [31] to the 
Schlesinger system of an arbitrary size. A similar structure of Hamiltonians can be found in a 'de- 
generate' version of the Gamier system studied by Kimura [33] and Shimomura [M]. Actually, 
this system coincides with the PI hierarchy associated with the string equation of type (2, 5). 

This paper is organized as follows. In Section 2, we introduce the string equations of type 
(2,p), and explain why they can be viewed as a higher order analogues of the first Painleve 
equations. In Section 3, these equations are cast into a 2 x 2 matrix Lax equation. Section 4 
is a brief review of the KdV and KP hierarchies. In Section 5, the PI hierarchy is formulated 
as a reduction of the KP (or KdV) hierarchy, and converted to 2 x 2 matrix Lax equations. In 
Section 6, we introduce the notion of spectral curve and consider the structure of its defining 
equation in detail. Though the results of this section appear to be rather technical, they are 
crucial to the description of Hamiltonians in spectral Darboux coordinates. Section 7 deals with 
the Hamiltonian structure of the Lax equations. In Section 8, we introduce spectral Darboux 
coordinates, and in Section 9, identify the Hamiltonians in these coordinates. In Section 10, 
these results are illustrated for the first three cases of (q,p) = (2, 3), (2, 5), (2, 7). 

2 String equation as higher order PI equation 

Let q and p be a pair of coprime positive integers. The string equation of type (q,p) takes the 
commutator form [4] 



K = H + AH. 



[Q,P] = l 
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for a pair of ordinary differential operators 

Q = dl+g 2 dl- 2 + ---+g q , P = dl + f 2 dl- 2 + --- + f p 

of order q and p in one-dimensional spatial variable x (d x = d/dx). In the following, we consider 
the equation of type (q,p) = (2, 2g + 1), g = 1, 2, . . .. We shall see that g is equal to the genus 
of an underlying algebraic curve (spectral curve). 

The simplest case, i.e., (q,p) = (2,3), consists of operators of the form 

3 3 

Q = dl + u, P = dl + -ud x + -u x , 
where the subscript means a derivative as 



du d 2 u 

OX OX z 



The string equation (|2.ip for these operators reduces to the third-order equation 
1 3 

jUxxx + -UU X + 1 = 0. 

We can integrate it once, eliminating the integration constant by shifting x — ► x + const, and 
obtain the first Painleve equation 

ju xx + -u 2 + x = 0. (2.2) 

The setup for the general case of type (2, 2g + 1) relies on the techniques originally developed 
for the KdV hierarchy and its generalization |35|, [36l [371 EH] • The basic tools are the fractional 
powers 

Q n+1/2 = + ^f^f- 1 + • • • + Rn+ldx- 1 + ■ ■ 

of Q = d 2 + u. The fractional powers are realized as pseudo-differential operators. The coeffi- 
cient R n +i of <9~ x is a differential polynomials of u called the Gelfand-Dickey polynomial: 

R = 1, R\ = — , R2 = -u xx + -u 2 , 
2 00 



1 3 



For all n's, the highest order term in R n is linear and proportional to u^ 2n 2 ^: 
1 



Rn 



,(2n-2) 



22n+2 

As in the construction of the KdV hierarchy, we introduce the differential operators 
B 2n+ i = (Q n+1/2 )> < n = 0,l,.... 

where ( )>o stands for the projection of a pseudo-differential operator to a differential operator. 
Similarly, we use the notation ( )<q for the projection to the part of negative powers of d x , i.e., 



X "J ''- =^2 a i 9 ^ W2 a 3 d l) =J2 a i dJ v 
\jez / >o i>o Vjez / <o j<o 
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These differential operators have the special property that the commutator with Q is of order 
zero. More precisely, we have the identity 

[P>2n+liQ] = %Rn+i,x- (2.3) 
Consequently, if we choose P to be a linear combination of these operators as 

P = B 2g+1 + Cl B 2g -i + ■■■ + c g B x (2.4) 
with constant coefficients c\, . . . , c g , the commutator with Q reads 

[Q, P] = -2R g+1)X - 2 Cl R g>x 2c g R ljX . 

The string equation (|2.ip thus reduces to 

2R g+hx + 2 Cl R g>x + • • • + 2c g R hx + 1 = 0. (2.5) 
This equation can be integrated to become the equation 

2R g+1 + 2 Cl R g + • • • + 2c g R 1 + x = 0, (2.6) 



which gives a higher order analogue (of order 2g) of the first Painleve equation 

When we consider a hierarchy of commuting flows that preserve the string equation, the 
coefficients ci,...,c g play the role of time variables. For the moment, they are treated as 
constants. 

3 Matrix Lax formalism of string equation 

The string equation (|2.ip is accompanied with the auxiliary linear problem 

= AV, Pip = d x iP (3.1) 

with a spectral parameter A (d\ = d/d\). We can rewrite this linear problem to a 2 x 2 matrix 
form [5] using various properties of the operators B>2 n +i and the Gelfand-Dickey polynomials R n 
[351 EH E3 EH] as follows. 

The first equation of (|3.1f) can be readily converted to the matrix form 

d x iP = C/ (A)V, (3.2) 

where 

To rewrite the second equation of (|3.ip . we use the 'Q-adic' expansion formula 

B>2n+1 = yRmdx — -R m ,x\Q n m (3.3) 

m=0 ^ ' 

of R>2n+iS- By this formula, i^n+iV' becomes a linear combination of ip and ip x as 

B>2n+\1p = ^2 [ Rm ^ x ~ 2 Rm ' x ^) = R m(X)^x ~ -Rm(X) x 1p, 
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where R n (X) stand for the auxiliary polynomials 

R n (X) = X n + R 1 X n ~ 1 + --- + R n , 71 = 0,1,..., (3.4) 

that play a central role throughout this paper. Since P is a linear combination of i^n+i's as (I2.4D 
shows, we can express Pip as 

PiP = a(A)V + /3(A)^x 
with the coefficients 

/3(A) = i? 9 (A)+c 1 ii 3 _ 1 (A) + --- + C9 i?o(A), a(A) = -i/3(A) x . (3.5) 

Moreover, differentiating this equation by x and substituting ip xx = (\—u)tp, we can express Pipx 
as 

P^ x = (P^) x = 7 (A)^ - a(X) x ^ x , 

where 

7 (A) = a(X) x + (A - u)/3(A) = -\f3{X) xx + (A - «)/3(A). (3.6) 

The second equation of the auxiliary linear problem (|3.ip can be thus converted to the matrix 
form 

d x iP = V(X)tp (3.7) 
with the coefficient matrix 



V(X) 



a(A) /3(A) 
7(A) -a(A) 



The differential equation (|3.2p in a: may be thought of as defining isomonodromic deforma- 
tions of the matrix ODE (|3 . Tf) in A with polynomial coefficients. The associated Lax equation 
reads 

d x V(\) = [U (\),V(\)] + U>(\), (3.8) 
where the last term stands for the A-derivative of Uq(X), i.e., 

U^X) = d x U (X) = (J J 

The presence of such an extra term other than a matrix commutator is a common characteristic 
of Lax equations for isomonodromic deformations in general. 

In components, the Lax equation consists of the three equations 

^a(A)= 7 (A)-(A-n)/3(A), d x /3(X) = -2a(A), ^ 7 (A) = 2(A - u)a(A) + 1. 

The first and second equations can be solved for a(A) and 7(A); the outcome is just the definition 
of these polynomials in (|3.5[) and fl3.6f) . Upon eliminating a(X) and 7(A) by these equations, the 
third equation turns into the equation 

-P(X) XXX - 2(A - u)/3(A) x + ^/3(A) + 1 = (3.9) 
for /3(A) only. 



Hamiltonian Structure of PI Hierarchy 



7 



Let us examine (j3.9|) in more detail. By the definition in (|3.5p . /5(A) is a polynomial of the 
form 

/3(A) = A 9 + Z^- 1 + • • • + p g 
with the coefficients 

Pn = Rn + c\R n -i H hc n i? , n=l,...,g. (3.10) 

Upon expanded into powers of A, (|3.9p yield the equations 

^(3 n ,xxx + 2uf3 n , x + u x /3 n - 2p n+ljX = 0, n = 1, . . . ,g - 1, 

and 

^P g ,xxx + 2u(3 9:X + u x (3 g + 1 = 0. 
Recalling here the well known Lenard recursion formula 

Rx,n+1 — ~~^Rn,xxx uR njX + ~U x R n , (^•H) 

one can see that the forgoing equations for n = 1, . . . ,g — 1 are identities, and that the latter 
one is equivalent to (|2.5D . 

4 KdV and KP hierarchies 

We shall derive the PI hierarchy from the KdV or KP hierarchy. Let us briefly review some 
basic stuff of these hierarchies [131 [T4l [T5] . Of particular importance is the notion of the Orlov- 
Schulman operator [16] that supplements the usual Lax formalism of these hierarchies. 

4.1 KdV hierarchy from KP hierarchy 

Let L denote the Lax operator 

L = d x + u 2 d~ 1 + u 3 d~ 2 H 

of the KP hierarchy. L obeys the Lax equations 

d n L = [B n ,L], B n = (L n )> , 

in an infinite number of time variables ti, £2, • • • (d n = d/dt n ). As usual, we identify t\ with x. 
The KdV hierarchy can be derived from the KP hierarchy by imposing the constraint 

(L 2 ) <0 = 0. 

Under this constraint, Q = L 2 becomes a differential operator of the form 

Q = d 2 + u, u = 2u2, 
all even flows are trivial in the sense that 

d2nL = [L 2n , L] = 0, 
and the whole hierarchy reduces to the Lax equations 

d 2n+1 Q = [B 2n+1 ,Q], B 2n+1 = (Q n+1/2 )> , (4.1) 
of the KdV hierarchy. By (|2.3|) , these Lax equations further reduce to the evolution equations 

d 2n+1 u = 2R n+ i :X (4.2) 

for u. 
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4.2 Orlov— Schulman operator 

The Orlov-Schulman operator is a pseudo-differential operator (of infinite order) of the form 

oo oo 

M = nt n L n ~ l +x + J2 VnL- n ~\ 

n=2 n=l 

that obeys the Lax equations 

d n M = [B n ,M] (4.3) 
and the commutation relation 

[L,M] = 1. (4.4) 

The existence of such an operator can be explained in the language of the auxiliary linear 
system 

Lip = zip, d n tp = B n ip, (4.5) 
of the KP hierarchy. This linear system has a (formal) solution of the form 

tp= ^l + ^WjZ-'j exp^xz + jr^nz 1 ^, (4.6) 

One can rewrite this solution as 



/ oo 



ip = W exp \xz + 2_^ t n z n , 

n=2 / 



where TV is a pseudo-differential operator (called the 'dressing operator' or the Sato-Wilson 
operator) of the form 



W = l + ^Wjd, 
i=i 



j 

X 



that satisfies the evolution equations 

d n W = -(W -d^-W-^^W (4.7) 
and the algebraic relations 

L = W ■ d x ■ W~\ B n = {W-3 n x - W- 1 )^. 
One can now define M as 

M = W {^ntnd™- 1 + W~ l , 

which turns out to satisfy the foregoing equations (|4.3p . (|4.4p and the auxiliary linear equation 
Mip = d z ip. (4.8) 

Remark 1. In the reduction to the KdV hierarchy, the 9~ 1 -part of (|4,7p gives the equation 
d 2n +iwi = -R n +i- (4.9) 

Since u = —2w\^ x , this equation may be thought of as a once-integrated form of (|4.2p . 
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5 PI hierarchy 

We now formulate the PI hierarchy as a reduction of the KP hierarchy, and rewrite it to a 2 x 2 
matrix Lax equation, as we have done for the string equation itself. The string equation and the 
commuting flows of the PI hierarchy are thus unified to a system of multi-time isomonodromic 
deformations. 

5.1 PI hierarchy from KP hierarchy 

To derive the string equation from the KP hierarchy, we impose the constraints [9] 

(Q) <0 =0, (P)<o = (5.1) 
on the operators 

1 oo oo 

Q = L 2 , P = -ML~ X = nt n L n ~ 2 + v n L- n ~ 2 . 

n=l n=l 

This leads to the following consequences. 



1. As it follows from the commutation relation (14, 4j) of L and M, these operators obey the 
commutation relation [Q,P] = 1. 

2. Under the first constraint of (|5.ip . Q becomes the Lax operator d 2 +u of the KdV hierarchy. 
In the following, as usual, we suppress the time variables with even indices, i.e., 

t 2 = u = ■ ■ ■ = o. 

3. The second constraint of (|5.ip implies that P is a differential operator (of infinite order) 
of the form 



1 °° n ,-. 



2 



Therefore, if we set 



2 

t2g+3 = 2g + \ '' t2g+5 = t2g+7 = ' ' ' = ' ^' 2 ' 

we are left with a differential operator of the form (12. 4p with the coefficients c\ , . . . , c g 
depending on the time variables as 

c n = c n (t) = 2n J" - frfo+i) n = l,...,g. (5.3) 



4. The auxiliary linear equations (|4.5p and (|4.8p imply the linear equations 

Qip = z 2 tP, = \z~ l d z *l>. 

These linear equations can be identified with the auxiliary linear equations (|3.ip of the 
string equation if we define A as 
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We can thus recover, from the KP hierarchy, the string equation (|2.ip of type (2, 2g + l) along 
with g extra commuting flows 

d2n+iQ = [B 2n +i,Q], d2n+iP = [B 2n+ i,P], n = l,...,g. 

We call this system the PI hierarchy. Let us mention that this hierarchy was first discovered in 
a more direct way [3J. Compared with that approach, the approach from the KP hierarchy [9] 
is more transparent. 

Remark 2. This hierarchy thus contains g + 1 flows with time variables t\(= x), £3, . . . , t 2g +i- 
The last flow in £23+1, however, turns out to be spurious, i.e., can be absorbed by other flows 
(see below). We shall suppress this flow when we consider the Hamiltonian structure. 

5.2 Matrix Lax formalism of commuting flows 

We can again use the Q-adic expansion formula (|3.3|) of B 2n +\ to rewrite the auxiliary linear 
equations 

d 2n +iip = B 2n +iip, n = l,...,g, 
of the PI hierarchy to linear equations 

d2n+iip = £4(A)V (5.4) 
for rj). The matrix elements of the coefficient matrix 
_ / A n (A) B n (A) 

Un[A) -{r n (\) — A n (A) 

are the following polynomials in A: 

B n (A) = Rn(X), A n (A) = ~Rn(X) x , 

T n (A) = ~Ra(X) xx + (A - u)R n (X). (5.5) 

(|3.2p can be included in these linear equations as a special case with n = (t\ = x). 

As one can see by comparing these matrix elements with those of C/ n (A)'s defined in ([3.5|) 
and (|3.6p . V(A) is a linear combination of U n (X)'s: 

V(X) = U g (X) + Cl (t)U g -x(X) + ■■■ + c g (t)Uo(X). (5.6) 

We shall consider implications of this linear relation later on. 

Having obtained the full set of auxiliary linear equations (|3.2p , (13. 7|) , (15. 4p in a 2 x 2 matrix 
form, we can now reformulate the commuting flows of the PI hierarchy as the isomonodromic 
matrix Lax equations 

d 2n+1 V(X) = [U n (X),V(X)] + U^X), n = 0,l,...,g, (5.7) 

including (13. 8p as the case with n = (t\ = x). Note that the Lax equations of the higher flows, 
too, have an extra term 

u' n (x) = d x u n (x) 

other than the matrix commutator [U n (X), V(A)]. 
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5.3 Variant of self-similarity 

Nowadays it is widely known that many isomonodromic equations can be derived from soliton 
equations as 'self-similar reduction' |39} I40j. The transition from the KP hierarchy to the PI 
hierarchy is actually a variant of self-similar reduction. 

Recall that P is a linear combination of B n 's as fj2.4f) shows. This implies that the left hand 
side of the auxiliary linear equation Pip = d\ip can be expressed as 

= (d 2g+1 + cx(t)d2g-i + ■■■ + c s (t)di)V>. 

The auxiliary linear equation thus turns into a kind of linear constraint: 

(d 2g+ i + ci(t)5 2fl _! + • • • + Cgit)^ = d x ^. (5.8) 

If we were considering the equation 

[Q,P}=0 

of commuting pair of differential operators [101 [HI H], all c n 's would be constants, and (|5.8p 
would mean the existence of stationary directions (in other words, translational symmetries) in 
the whole time evolutions of the KdV hierarchy. It is well known that this condition characterizes 
algebro-geometric solutions of the KdV hierarchy. 

In the present setting, where c n 's are not constant but variables as (|5.3p shows, (|5.8p may be 
thought of as a variant of self-similarity condition. We can see some more manifestation of this 
condition. For instance, V(A) satisfies the linear equation 

(d 2g+l + Cl {t)d 2g -i + ■■■ + Cg (t)9i)y(A) = v'(\) 

as one can deduce from the Lax equations (|5.7p and the linear relation (15. 6|) among their 
coefficients. Similarly, u satisfies a similar equation 

{d 2g+1 + ci{t)d 2g -i H h c g {t)di)u = 1 

as a consequence of (I2.5P and (|4.2p . These equations show, in particular, that the t 5 +i-flow is 
spurious. 

Remark 3. (|5.8p stems from the Virasoro constraints associated with the string equation. The 
Virasoro constraints are usually formulated in the language of the r-function [6] . Reformulated in 
terms of ip, the lowest one (the L_i-constraint) of those constraints becomes the linear equation 

^ — ^— t 2n+ id 2n -itp = d\ip. 

n=0 

This turns into the aforementioned constraint (|5.8p if the higher time variables £23+3 ; ^2g+5 > • • • 
are set to the special values of (|5.2p . 

6 Building blocks of spectral curve 
6.1 Equation of spectral curve 

We now consider the spectral curve of the matrix V(A). This curve is defined by the characte- 
ristic equation 

det(//J - V(A)) = n 2 + det V{\) = 
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or, more explicitly, by the equation 

/j 2 = h(X) = a(A) 2 + 0(X)j(X). 
Since q(A), /3(A), 7(A) are polynomials of the form 

a(A) = aiA 9_1 + • • • + a g , 

/3(A) = \9 + + ■ ■ ■ + (3 g , 

7 (A) =A^ 1 + 7l A3 + ... + 75+1 , 

h(X) is a polynomial of the form h(X) = X 2g+1 + ■ • • , and the spectral curve is a hyperelliptic 
curve of genus g. 

The spectral curve plays a central role in the algebro-geometric theory of commuting pairs 
I12j. In that case, the flows of the KdV hierarchy can be translated to the isospectral Lax 
equations 

d 2n+1 V(X) = [U n (X),V(X)] (6.1) 

of the same matrix V^(A) as we have used thus far (except that c n 's are genuine constant). 
In particular, the polynomial h(X) (hence the spectral curve itself) is invariant under time 
evolutions: 

d 2n+1 h(X) = 0. 

This system (|6.ip is called 'the Mumford system' |17j . 

In contrast, the polynomial h(X) for the matrix Lax equations (15, 7ft of the PI hierarchy is 
not constant in time evolutions. Straightforward calculations show that the t-derivatives of h(X) 
take non-zero values as 

d 2n+l h(X) = TiU' n (X)V(X). 

For instance, the lowest equation for n = (t± = x) reads 

d x h(X) = (3(X). 

The spectral curve thus deforms as x and i2n+i's vary. Obviously, the extra terms U' n {X) are 
responsible for this phenomena. This is a common feature of isomonodromic Lax equations. 



6.2 Some technical lemmas 

As we shall show later on, the coefficients of terms of higher degrees in the polynomial h{X) 
have a rather special structure. We present here a few technical lemmas that are used to explain 
this fact. 

The following lemmas are concerned with the KdV hierarchy rather than the PI hierarchy. 
In the case of the KdV hierarchy, the matrices U n (X) are defined for all nonnegative integers 
n = 0,l,2,.... 



$(A) 



Lemma 1. There is a 2 x 2 matrix 

1 + OiX- 1 ) oix- 1 ) 
w^OiX" 1 ) 1 + 0(X~ 1 ) 

of Laurent series of X that satisfies the equations 

d 2n+1 <S>(X) = C/„(A)$(A) - $(A)A n A, n = 0,l,2,..., (6.2) 

where 



A 



1 

A 
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Proof. Let ip{z) be the special solution, (|4.6p . of the auxiliary linear equations. In the case of 
the KdV hierarchy, this is a function of the form if)(z) = w(z)e^ z \ where 

oo oo 

w(z) = 1 + Y^ Wj z-i, £{z) =Y, t 2n+iz 2n+1 (h = X). 

j=l n=0 

The associated vector-valued function 

m = ( f/f ) = ( , f Z) < ^ ) ^ 



jl>{z) x ) \zw(z) + w(z) 

satisfies the auxiliary linear equations 

d 2 n+Mz) = U n (X)ip(z) (A = z 2 ) 

with the same coefficients U n (X) as in (|5.4[) but now defined by (|5.5|) for all nonnegative integers 
n = 0,l,2,.... 

Moreover, since A remains invariant by substituting z — > —z, i/)(—z) is also a solution of these 
linear equations. Thus we actually have a 2 x 2 matrix- valued solution 

(it \ it- \\— ( w ^ w (~ z ) \ ( ei{z) 

{V{Z) VK Z) >~ \zw{z)+w{z) x -zw{-z)+w{-z) x ) V e-«*) 



We now consider 

¥(A) = (i>(z) i/>(-z)) 



1 1\ 1 ( ip(z) +ip(-z) ip(z) -if>{-z) 



z -z \ 2 2z 



which is also a matrix- valued solution of the foregoing linear equations. Note that the matrix 
elements are even functions of z (hence functions of A). Moreover, ^(A) can be factorized to the 
product of 

= ( w{z) w{ ~ z) V 1 l Y 

\zw(z)+w(z) x —zw(—z)+w(—z) x J \z —zj 

and 

C-')c: ^K-r-MP^ + ) 

By plugging 

*(A) = <&(A)exp(f>„ +1 A 2 " +1 J 

\n=0 / 

into the auxiliary linear equations, $(A) turns out to satisfy the equations 

d 2n+1 $(X) = C/- n (A)*(A) - $(A)A 2 " +1 , 

which are nothing but (|6.2j) because of the identity A 2n+1 = X n A. ■ 
Now let us introduce the matrix 

U(X) = $(A)A$(A) _1 . 
As a consequence of (|6.2p . it satisfies the Lax equations 

d 2n+1 U(X) = [U n (X),U(X)}, n = 0,1,2,.... (6.3) 
The following lemma shows that we can use U(X) as a kind of generating function for U n (Xys 
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Lemma 2. The matrix elements of 

U{X) -{r(X) -A(A)J' 
are Laurent series of X of the form 

A(A) = 0(A^), B(A) = l + OiX" 1 ), r(A) = A + 0(A°), 
that satisfy the following algebraic conditions: 

A n (X) = (A"A(A))> , B„(A) = (A»B(A))> , 

r n (A) = {X r T(X))> -R n+1 , n = 0,l,..., (6.4) 
A(A) 2 + B(A)r(A) = A, (6.5) 

where ( )>o stands for the polynomial part of a Laurent series: 

J£Z ) > 3>0 



Proof. Rewrite (|6.2p as 

U n (X) = d 2n+l $(X) ■ ^(Xy 1 + cD(A)A n A<I>(A)- 1 = d 2n+1 3?(A) ■ ^A)" 1 + A n *7(A), 

and compare the polynomial part of both hand sides. U n (X) is a matrix of polynomials, and the 
first term on the right hand side is a matrix of the form 



d 2n+ iH\) ■ Hxy 1 



0{X^) 0(X~ 1 

d 2n+lWl + o(A- 1 ) o(x~ l ) 



Thus, also by recalling (|4.9p . (|6.4p turns out to hold. (|6.5|) is an immediate consequence of the 
definition of U(X) and the identity det A = A. ■ 

Remark 4. Since B n (A) is equal to the auxiliary polynomial R n (X) defined in (|3.4p . the second 
equation of (|6.4p implies that B(A) is a generating function of all i? n 's: 

B(A) = l + J RiA" 1 + J R 2 A~ 2 + -- - . 

The lowest (n = 0) part of (|6.3|) is not a genuine evolution equation. In components, this 
equation reads 

d x A(X) = r(A) - (A - u)B(X), 8 X B(X) = -2A(A), 

d x T(X) = 2{X-u)A(X). (6.6) 
The first two equations can be solved for A(A) and T(A) as 

A(A) = ~B(X) X , T(A) = ~B(X) XX + (A - u)B(A). 
The third equation thereby reduces to 

^B(X) XX - 2(A - u)B(A) x + n,B(A) = 0. (6.7) 
It is easy to see that this is a generating functional form of the Lenard relations (|3.1ip . 
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Remark 5. If A(A) and T(A) are eliminated by (|6.6p . (|6.5p becomes another generating func- 
tional formula 

\{B(X) X ) 2 + B(A) (~B(X) XX + (A - u)B(A)) = A (6.8) 

of relations among i? n 's. (|6.8[) may be thought of as a once-integrated form of (|6.7p . If expanded 
in powers of A, (16, 8h becomes a sequence of relations that determine R n recursively without 
integration procedure. This is an alternative and more practical way for calculating i? n 's. 

6.3 Detailed structure of h(X) 

We now turn to the issue of h(X). Let us recall (|5.6p . In components, it reads 

q(A) = A g (X) + ci(i)A a _i(A) + • • • + c g (t)A {X), 
/3(A) = B fl (A) + ci(t)B a _i(A) + • • • + c s (t)B (A), 

7 (A) = r 9 (A) + ci(t)r ff _!(A) + • • • + c 9 (t)r (A). 

We can thereby express /i(A) as 

9 

h(X)= c m (t)c n (t){A g _ m (X)A g _ n {X) + B 9 _ m (A)r 9 _ n (A)), 

m,n=0 

where it is understood that co(t) = 1. Moreover, since (|6.4p implies that 

A™ (A) = A n A(A) + ©(A" 1 ), B„(A) = A n B(A) + ©(A" 1 ), 
T n (A) =A n r(A)-i? n+1 + 0(A- 1 ), 

we can further rewrite h(X) as 

fc(A)= £ c m (t)c n (t)A 2 ^ m - n (A(A) 2 + B(A)r(A)) 

m,n=0 
9 

- 2 c m (t)R g+1 „ m X 9 + 0(A^ X ). 

m=0 

We can now use (|6.5p and (|2.6|) on the right hand side. This leads to the following final result. 
Theorem 1. Up to terms of 0(X 9 ~ 1 ), h(X) can be expressed as 

h(X) = X 29+l + 2 Cl (t)X 29 + (2c 2 (t) + ait^X 29 - 1 + (2ca(t) + 2 Cl (t)c 2 (i))A 29 - 2 + • • • 

+ ^ c m (t)c g . m (t)X 9+1 c m {t)c g+l _ m {t)+x\x 9 + 0{X 9 ~ 1 ). (6.9) 

m=0 \m=l / 

In particular, the coefficients of A 2fil+1 , . . . , X 9 in h(X) do not contain u, u x , .... 

Let Io(X) denote the part of h(X) consisting of A 29+1 , . . . , A 9 , and Ji, . . . , I g the coefficients 
of A 9 " 1 , ... ,1: 

h(X) =I (X)+I 1 X 9 ~ 1 + ■■■ + !,. 
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We have seen above that Io(X) is a kinematical quantity that is independent of the solution 
of the PI hierarchy in question. In contrast, the remaining coefficients I±, . . . ,I g are genuine 
dynamical quantities. 

In the case of the Mumford system (|6.1[) . these coefficients I±,...,I g are Hamiltonians of 
commuting flows. More precisely, it is not these coefficients but their suitable linear combina- 
tions Hi, . . . , H g that exactly correspond to the flows in ti, t3, . . .. We shall encounter the same 
problem in the case of the PI hierarchy. 



7 Hamiltonian structure of Lax equations 

We use the same Poisson structure as used for the Mumford system [18} 119] , This Poisson 
structure is defined on the 3g+ 1-dimensional moduli space of the matrix V(A) with coordinates 
ai, . . . , a g , . . . , 7i, . . . , 7 9 +i- It is customary to use the tensor notation 

{V(A) f Vfa)} = Yl { V <*W, V cd {\)}E ab ® E cd , 

a,b,c,d=l,2 

where E ab denote the usual basis of 2 x 2 matrices. The Poisson brackets of the matrix elements 
of V(A) can be thereby written in a compact form as 

{V(\) ® V(ji)} = [V(X) ® I+I V(ji),r(\-n))] + \V(\) ® /-/ ® V(p), E 2 i ® E 2 i], 

where r(X — fi) is the standard rational r-matrix 

p 

r(\-n) = j—, P= Yl E ab ®E ba . 

" a,6=l,2 

This is a version of the 'generalized linear brackets' [20]. More explicitly, 
WA),4)} = 0, {0(\),p(n)} = O, 

{a(A), (5 (/i)} = — , {a(A), 7 (/i } = , 

A — p A — /i 

{/?(A), 7 (/x)} = 2 a(A j~ a(/x) , { 7 (A), 7 (M)} = -2a(A) + 2a( M ). (7.1) 

We can convert the Lax equations (|5.7p of the PI hierarchy to a Hamiltonian form with 
respect to this Poisson structure. This procedure is fully parallel to the case of the Mumford 
system. 

A clue is the Poisson commutation relation 



{V(X),h(p)} 



V(X),^-+P(fi)E 2 i 
A — fx 



(7.2) 



which can be derived from (|7.1f) by straightforward calculations. One can derive from this 
relation the Poisson brackets {V(A), I n +i} a s follows. 

Lemma 3. 

{V(X),I n+ i} = [V n (X),V(X)], (7.3) 
where V^(A) is a matrix of the form 



V n (X) 



a n (X) p n (X) 
7 n (A) -a„(A) 



Hamiltonian Structure of PI Hierarchy 



17 



with the matrix elements 

a n (X) = (A"-%(A))> , (3 n (X) = {X n -^(X))> , 

7n(A) = (A"- 9 7 (A))> -/3 n+ i, n = 0,1,..., 3-1. 
Proof. Tn+i can be extracted from h(fj.) by a contour integral of the form 

I n+ i = f^ n - 9 h(»), 

where the contour is understood to be a circle around fi = oo. The same contour integral applied 
to (|7,2p yields the Poisson bracket in question: 

{V(X),I n+1 } = j> ^-9 {V{X)M ^ )} = [vi^f^ ^M +(3 n+ iE 21 

Let us examine the matrix 

J 2m A — fi 

Since this type of contour integral gives, up to signature, the polynomial part of a Laurent 
series /(A) as 

£££ = -(«*»*■ 

V^(A) can be expressed as 

V n (X) = {X n ~9V(X))> - (3 n+ iE 21 . 

The statement of the lemma follows from this expression of V n (X). ■ 
Actually, this result is not what we really want - we have to derive [U n (X), V(A)] rather than 
[Ki(A), V(A)]. Here we need another clue, which is the following linear relations among {7 n (A)'s 
and y n (A)'s. 

Lemma 4. 

V (X) = U (X), 

V n (X) = U n (X) + ci(t)U n -x(X) H h c n (t)U (X), n = l,...,g-l. (7.4) 

Proof. The linear relations f|3. 10|) among /3 n 's and -R n 's imply the linear relations 

n {X) = R n (X) + atyRn-iiX) +•■■ + C n (t)R (X) 

= B„(A) + c 1 (i)B n _ 1 (A) + • • • + c n (t)B (A) 
of the auxiliary polynomials. On the other hand, a(A) and 7(A) s are connected with /3(A) as 

a(A) = ~(3(X)x, 7(A) = -\(3{X) XX + (A - u)(3(X), 

Expanding these relations in powers of A and picking out the terms contained a n (X) and 7n(A), 
one can see that a n (X) and 7 n (A) are linearly related to A n (A)'s and T n (A)'s with the same 
coefficients as 

a n {X) = A„(A) + ci(i)A„_i(A) + • • • + c n (t)A (A), 

7n(A) = r n (A) + ci(t)r re _i(A) + • • • + c n (t)T (X). 

These are exactly the linear relations presented in (|7.4p in a matrix form. ■ 
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In view of this lemma, we define new Hamiltonians H\, . . . ,H g by the (triangular) linear 
equations 

h = H h I n+1 = H n+1 + c 1 (t)H n + --- + c n (t)H 1 , n=l,...,g-l. (7.5) 

Note that H g+ i is not defined (because I g+ \ does not exist). The foregoing formula (|7.3p of the 
Poisson brackets of V(X) and I n 's can be thereby converted to the form 

{V(X),H n+1 } = [U n (X),V(X)} 

that we have sought for. We thus eventually obtain the following result. 

Theorem 2. Except for the t2 g +\-flow, the matrix Lax equations (|5.7p of the PI hierarchy can 
be cast into the Hamiltonian form 

d 2n+l V{X) = {V{X),H n+l } + U' n {X) 1 n = 0,1,... (7.6) 

with the Hamiltonians defined by (|7.5f) . 

Remark 6. As regards the excluded t2g+i-fl° w i the polynomial h(X) obviously contains no 
candidate of Hamiltonian. If we naively extrapolate (|7.5p to n = g, we end up with the linear 
relation 

= Ig + l = H g + 1 + Cl{t)Hg + ■■■ + Cg^Hi 

of the Hamiltonians. In a sense, this is a correct statement, which says that H g+ \ is not an 
independent Hamiltonian. 

Remark 7. Io(X) is a central element (i.e., a Casimir function) of the Poisson algebra. To see 
this, note that the right hand side of (|7.2p is of order 0(fi 9 ~ 1 ) as \i — > co. This implies that the 
terms of degree greater than g in h(fi) have no contribution to {V (X) , h(/i)} , in other words, 

{V(X),I o (ji)} = 0. (7.7) 

This means that Io(fi) is a Casimir function. This fact is in accord with the observation in the 
last section that Io(X) does not contain genuine dynamical variables. 

Remark 8. I n 's are Poisson-commuting, i.e., {Ij,Ik} = 0. This is a consequence of another 
basic Poisson relation 

{h(X),h(f,)} = 0, 

which, too, can be derived from (|7.ip by straightforward calculations (or by a standard r-matrix 
technique) . 

8 Spectral Darboux coordinates 

The construction of 'Spectral Darboux coordinates' is also parallel to the case of the Mumford 
system. These coordinates consist of the roots Ai, . . . , X g of /3(A) and the values /ii, . . . , fx g of 
a(A) at these roots of /3(A): 

9 

/3(A) = Y[(X - Xj), fij = a(Xj), j = l,...,g. 
i=i 

To avoid delicate problems, the following consideration is limited to a domain of the phase space 
where Aj's are distinct. 
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Xj and /ij satisfy the equation of the spectral curve: 



We thus have a g-tuple (Xj, fJ<j)j =1 of points of the spectral curve (in other words, an effective 
divisor of degree g) that represents a point of the Jacobi variety of the spectral curve. In 
the case of the Mumford system, the commuting flows are thereby mapped to linear flows on 
the Jacobi variety [111 W2\ [T7] . The case of the PI hierarchy is more complicated because the 
spectral curve itself is dynamical. If one wishes to pursue this approach, one has to consider 
the coupled dynamics of both the divisor and the underlying spectral curve; unlike the case of 
isospectral problems, this does not reduce the complexity of dynamics of the original nonlinear 
problem. Actually, this is not what we seek for. We simply borrow the idea of spectral Darboux 
coordinates to describe the Hamiltonian structure of the system in question. 

As it follows from (|7. If) by a standard procedure [20], [261 EZ] , these new variables satisfy the 
canonical Poisson relations 

{Aj,A fc } = 0, {fJ>j,Hk} = 0, {Xj^k} = $jk- 

On the other hand, they Poisson-commute with Iq(X), 

{A i ,/ (A)} = 0, {/i i ,/ (A)} = 0, 

because Io(X) is a Casimir function as (|7.7p shows. Thus Aj's and ^'s may be literally called 
'Darboux coordinates'. 

These Darboux coordinates Xj, [ij and the coefficients of X 2g , . . . , X 9 in Io(X) = X 2g+1 + • • • 
give an alternative (local) coordinate system of the 3g + 1-dimensional Poisson structure on the 
space of L-matrices with the original g + g + (g + 1) coordinates jj, (3j, jj. To reconstruct the 
L-matrix V(X) from these new coordinates, we use the familiar Lagrange interpolation formula 



3- 



^>{Xj)X-Xj 

that holds for any polynomial f(X) = f\X 9 ~ l + ■ ■ ■ + f g of degree less than g. Since 



/3(A) d(3(X) 



X — Xn 



dX, 



Edfin . g-. n 
9A~7 ' 

n=l J 



this formula implies the formula 



t /9'(Ai) dXj 



for the coefficients of /(A) as well. Note that f3 n 's are understood here to be functions of Aj's 
(in fact, they are elementary symmetric functions). We apply this formula (|8.2p to a(X) and 
obtain the explicit formula 



a. 



Y- 



fri P'(Xj) dXj 

that recovers a n 's from Aj's and /Ltj's. In a similar way, we apply ([8.2[) to the case where 



/(A) = ^I„A^ = MA)-/o(A), 



71=1 
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and find the expression 

J --^-^rOA- (8 ' 3) 

of I n 's in terms of Ay's, fj,/s and io(A). For instance, 
_ « /x 2 - IpjXj) 

In these formulas, Io(X) is understood to be the polynomial 

J (A) = A 29+1 + 2 Cl (t)X 29 + (2c 2 (i) + ci(t) 2 )A 21 '- 1 + (2c 3 (t) + 2 Cl (i)c 2 (t))A 25 - 2 + • • ■ 

9 /9~1 \ 

+ C m (t)Cg- m (t)\ 9 + 1 + C m (t)Cg+l_ m (t) + X \X 9 . 

m=0 \m=l / 

Once ct(A) and /i(A) are thus reconstructed, we can recover 7(A) as 
h(X) -a(X) 2 



7(A) 



/3(A) 



It is convenient to rewrite the foregoing formula ()8.3[) slightly. Recall the auxiliary polyno- 
mials 

f3 n (X) = X n + f3 1 X n ~ 1 + ■•• + /?„. 
Lemma 5. 



5A^ 



&_i(Aj), n = !,...,<?. (8.4) 



Proof. Start from the identity 

em = m = i( A ) - ^( A i) 

A — Xj A — Aj 

and do substitution 



-Pn-, r 1 - = -Pn(X + A.-A"" 2 + • • • + A"- 2 A + A"" 1 ) 

A - A j 



for each term on the right hand side. This leads to the identity 
00(A) 



dX 



j 



-- -X a ~ x - (X, } + /3!)A^ 2 (Af 1 + ftAf 2 + ■ ■ ■ + /3 fl _i) 

. A9 -i_ A(A . )A9 -2 ^(Aj), 



which implies (|8.4p . 

By these identities, we can rewrite (18.31) as 

^n+l=V ^77T-T Pn(Aj), « = 0, 1, . . . , g - 1. 
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(For notational convenience, n is shifted by one). We can derive, from these formulas, a similar 
expression of the Hamiltonians H n+ \ introduced in the last section. Recall that /3 n 's are con- 
nected with i? n 's by the linear relation (|3.10p . It is easy to see that the auxiliary polynomials 
(3 n (X), too, are linearly related to the auxiliary polynomials R n (X) as 

(3 n (\) = R n (X) + Cl(i)i?n_l(A) + • ■ ■ + C n (t)R (X). 

Comparing this linear relation with the linear relation (17. 5p among I n 's and i/ n 's, we find that 
the Hamiltonians H n+ \ can be expressed as 

H n +i = > ai( \\ — R n( X j)' n = 0, 

Note that i? n 's in this formula have to be redefined as functions of Xj that satisfy the linear 
relations (|3. 1U|) : 

Rl=Pl-Cl(t), R 2 =P2 -CiWfa +C l (tf -C 2 (t), .... 

In particular, these Hamiltonians are time-dependent, the time-dependence stemming from both 
J (A) and Rn(X). 

This is, however, not the end of the story. As it turns out below, these Hamiltonians (except 
H\ = Ji) do not give correct equations of motion in the Darboux coordinates Xj, fij. Correct 
Hamiltonians are obtained by adding correction terms to H^s. 

9 Hamiltonians in Darboux coordinates 
9.1 Equations of motion in Darboux coordinates 

Let us derive equations of motion for Xj's and //j's from the Lax equations (|5.7p . In components, 
the Lax equations take the following form: 

d 2n+l a(X) = B n (A) 7 (A) - P(X)T n (X) + A' n (A), 
d 2 n+i(3(X) = 2A n (A)/3(A) - 2B n (A)a(A) +B^(A), 

d 2n+ MX) = 2T n (X)a(X) - 2A n (A) 7 (A) + r' n (A). (9.1) 

To derive equations of motion for Ay's, we differentiate the identity f3{Xj) = by t 2n+ \. By 
the chain rule, this yields the equation 

d 2n +iP(\)\\=\ j +p'(x J )d 2n+1 x J = o. 

By the second equation of (|9.ip . the first term on the right hand side can be expressed as 

WWUj = -2B n (A)a(A j ) + B^(Aj) = -2^B n (A) + B' n (Xj). 
Thus the following equations are obtained for Ay's: 

d ^ x >- pfi) ~y(x~j- (9 - 2) 

To derive equations of motion for fJ,j's, we differentiate fj,j = a(Xj) by t 2n+ \. The outcome is 
the equation 

d 2n +ifJ-j = d 2n+ ia(X)\x=\ 3 + a'(Xj)d2n+\Xj. 



22 



K. Takasaki 



By the first equation of (|9.ip . the first term on the right hand side can be expressed as 
d 2 n + ia(X)\ x =x ] = B n (Aj)7(Aj) + A' n (Xj). 

The derivative d 2n +iXj in the second term can be eliminated by (|9.2p . We can thus rewrite the 
foregoing equation as 



2 H gi{\ j )+{3\\ ] ) 1 {\ ] ) a' (Xj)B' n (Xj) s 

&2n+l ^ ~ WW i j) WW) n{ >' 



Note here that the numerator of the first term on the right hand side is just the value of h!{\) 
at A = Aj: 

h'(Xj) = 2a(X J )a'(X j ) + /3'(A J ) 7 (A J ) + /?(A J ) 7 , (A J ) = 2 Mi «'(A i ) + ^'(A i ) 7 (A i ). 
Thus the following equations of motion are obtained for fi^s: 
_ h'(Xj)B n (Xj) a' (Xj)B' n (Xj) 

d2n+i ^-^w www - +K{X]) - (9 - 3) 

9.2 Why Hamiltonians need corrections 

We now start from the Hamiltonian form (|7.6p of the Lax equations and repeat similar calcu- 
lations. In the case of isospectral Lax equations, such as the Mumford system, this procedure 
should lead to a Hamiltonian form of equations of motion for the spectral Darboux coordinates. 
In components, (|7.6p consist of the following three sets of equations: 

cWWA) = {a(A), H n+1 } + A'„(A), 
<Wi/?(A) = {/3(A), H n+1 } + B' n (A), 

d2n+il(X) = { 7 (A), H n+1 } + r;(A). (9.4) 
We again start from the identity 

= d 2n+1 f3(X j ) = d 2n+ il3(X)\x=x J +P'(X j )d 2n+l X j 
and consider the second equation of (|9.4p . which implies that 

d 2n+ il3(X)\\=\ j = {/3(A),F„ + i}|a=a j +KW- 
To calculate the first term on the right hand side, we use the identity 

= {(3{Xj), H n+ i} = {/3(A),# n+ i}| A=A , +p'(X j ){X j ,H n+1 }. 
Thus the following equations of motion are obtained for Aj's: 

d 2n+1 X, = {X„H n+1 }-^^ (9.5) 

In much the same way, we can derive the following equations of motion for ^j's: 

d 2n+m = fo, H n+1 } - + KW (9-6) 
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These results clearly show that H n+ \ is not a correct Hamiltonian. If we could find a correct 
Hamiltonian, say K n+ i, the equations of motion would take the canonical form 

C^n+lAj = {Xj,K n+ i}, (?2n+lMj = {Hj,Kn+l}- 



(|9.5p and (|9.6p fail to take this form because of extra terms on the right hand side. These terms 
stem from the extra term U' n {\) in the Lax equations (|5.7p . When we say (|5,7p is 'Hamiltonian', 
we ignore the presence of this term. It, however, cannot be ignored if we attempt to formulate 
the Hamiltonian structure in the language of the spectral Darboux coordinates Xj and [ij. 

The case of n = is exceptional. Since Bo (A) = 1, the extra terms on the right hand side 
of (19.21) and (19. 3D disappear. Therefore 



/ (Aj) 

is a correct Hamiltonian for the equations of motion for the ii-flow (namely, the higher order 
PI equation itself). Note, however, that the extra term Uq(X) in the Lax equation still persist. 

Remark 9. Another consequence of the foregoing calculations is that the Poisson brackets 
of H n+ i and the Darboux coordinates (in other words, the components of the Hamiltonian 
vector field of H n+ \) are given by 

2 H B n (X 3 ) fe / (A J -)B n (Aj) 
{Aj ' if " +l} " /3'(A,) ' " /3'(Aj) ■ 

Consequently, the Poisson brackets of I n +i and the Darboux coordinates turn out to be given by 

n r i _ 2^/3n(A J ) _ h'(Xj)(3 n (\j) 

I n+1{ — r,/f\ \ — ' lA t i) J n+l/ 



/?'(A,) ' /3'(A,) • 

One can derive these results directly from the Poisson brackets (|7.ip of the matrix elements 
of V(A) as well. 

9.3 Corrected Hamiltonians 

We now seek for correct Hamiltonians in such a form as 

Kn+i = H n+ \ + AH n+ \. 
Correction terms AH n+ i have to be chosen to satisfy the conditions 

{Xj,AH n+1 ^ - - p , {fj,j,AH n+ i} pfjyj h A n(Ajj- 

We can convert this problem to that of 7 n +i's, namely, the problem to identify the correction 
term to I n+ i's: 

AI n+1 = AH n+1 + Cl {t)AH n + ■■■ + Cn-^AHx. (9.7) 
The conditions for to A/ n+1 read 

r\ at \ ^n(Aj) r AT i <x'{^j)Pn( X j) , _/ /x N 

{Xj,Al n+1 \ - - p , {fj,j, AI n+1 \ ^(\~) + a ^ X i> 
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or, equivalently, 

dAI n+1 _ [3' n (\ 3 ] 



dAI n 



+1 



dfii 



0>{Xj 



ax* 



a' n (Xj). 



(9.8) 



This slightly simplifies the nature of the problem. 

The goal of the subsequent consideration is to prove that a correct answer to this question is 
given by 



n+1 



R'l 



k=l 



(9.9) 



Obviously, the first half of (19. 8ft is satisfied; what remains is to check the second half. 
We can use the following lemma to reduce the problem to each term of the sum in (|9.9 



Lemma 6. 



Y- 

Rl 



{P'{x k ) dx k 



fi k d(3' n (X) 



A=A, 



dX k 



A=A ; 



Proof. By the Lagrange interpolation formula (|8.ip . a(X) can be expressed as 



a(X) 



Y- 



fi k d(3(X) 



to both hand sides yields another identity 



k ^P'(X k ) dX k ■ 

Applying the projection operator (X n ~ 9 

r\\ d/3 n (X) 
^ (3' n (X k ) dX k 

The statement of the lemma follows by differentiating both hand sides of these identities and 
by setting A = Xj. ■ 
Because of this lemma, checking the second half of (|9.8|) can be reduced to proving the 
following identity: 



1 



dXj 0>{x k ) p(\ k ) dx k 



(3' n {\j) dp'(X) d(3' n {X) 



dX k 



(9.10) 



\=\j 



Note that this is a genuinely algebraic problem (related to elementary symmetric functions). 
We prepare some technical lemmas. 



Lemma 7. 

d 2 (5 n 



1 



dXjdX k 



Ai — At V 9 X 



d/3 n dp n 



dX k 



(j ¥= k), 



9*1 



0. 



(9.11) 



Proof. Differentiate 



dX k 



II( A - A 



once again by Xj. If j / k, this yields the identity 



dXjdX k 



I(A-A, 



(A — Xj)(X — Afe) 



1 fd(3(X) dp(X)\ 
Xj - A fc V dXj dX k J ' 



proving the first part of (|9.11|) . Similarly, if j = k, the outcome is the identity d 2 j3{X) / dX 2 , = 0, 
which implies the rest of (|9.11j) . ■ 
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Lemma 8. 

^ {Xk) --w k -^r Xk d\ k Xk ■ (9 - 12) 



Proof. (3' n (X) can be expressed as 

(3' n {\) = nX 1 - 1 + (n - l)/3iA™- 2 + • • • + /3 n _i 

= (A n_1 + /?iA n ~ 2 + • • • + /3 n _!) + • • • + (A + /?i)A n " 2 + A™" 1 
= /3 n _i(A)+/3 n _ 2 (A)A + --- + /? (A). 

Upon substituting A = X k and recalling (|8.4|) . (|9.12p follows. 
Lemma 9. 



a^(A fc ) i 



/^'(A fc ). (9.13) 



9A fe 2 

Proof. If j ^ k, (I9TT21) and f[97TT]) imply that 

-Au — ■ ■ ■ — — — — — At, 



dXj dXjdXk dXjdXk dXjdXt 



4 fe 



1 (d(3 n d(3 n \ + /fi_9ft\ r i 

\j-Xk\a\j dXk) \dXj ax k J k 

Xj — Afc V ^ A j 9Aj fc 5Aj fe 



1 



/ dj3n + A + h ^n-i 

A, - A fc V^^fc <9A fc 9A fc fc 



Obviously, 



<9/? n 9/3 n _i _! <9/3 n (A A 

H ^ — -U + • • • + TT"; Ai 



9Aj ' 9Aj " K ' ' dXj' k dXj ' 

and by (gmj) , 

9AT + ^AT Afe+ --- + ^ Afe =-^y 

Thus the first part of (p7[3j) follows. If j = k, (I9TT2D and (|9TTT|) imply that 

On the other hand, differentiating the identity 

^(A) = /3 n _!(A) + /? n _ 2 (A)A + • • • + /? (A)A n ~ 1 

(which has been used in the proof of (|9.12p ) yields 

^(A) = /3^_i(A) + &_ 2 (A)A + • • • + (5' 1 (X)X n ~ 1 

+ /3 n _ 2 (A) + 2/3 n _ 3 (A)A + ... + («- l)/?o(A)A™- 2 . 
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One can eliminate the derivatives (3' n _ 1 (X), . . . , (3[(X) by the preceding identity itself. The out- 
come reads 

P' r [(\) = 2(/3 n _ 2 (A) + 2/3 n _ 3 (A)A + • • • + (n - 1)/3 (A)A"- 2 ). 
Upon substituting A = X k and using (|8.4p . one finds that 

&'(A fc ) = 2(/3 n _ 2 (A fc ) + 2/? n _ 3 (A fc )A fc + • • • + (n - l)/? (A fc )A£- 2 ) 
= -2| 



^- 1 +2^A fc + 



<9A A 



dA* 



+ (n _ 1)1^-2 



The last identity and the foregoing expression of df3' n {X k ) / 'dX k lead to the second part 

of flajjjp . ■ 

Using these lemmas, we can calculate both hand sides of (|9.10p . 

Let us first consider the case of j 7^ k. The left hand side of (|9.10j) can be calculated by the 
Leibniz rule and (|9.13p . Note here that the formula (|9.13p for n = g takes such a form as 



dX j 



Xi — Xi 



because j3 g {X k ) = !3{X k ) = 0. The outcome of this calculation reads 

d PUXk) _ 1 d(5 n {X k ) 

dXj(3'(X k ) /3'(Xk)(Xj — Xk) dXj ' 

As regards the right hand side of (|9,10p . we can use (|9.13p to calculate the derivatives in the 
parentheses as 



d[3'(X) 



dX k 
d{5' n {X) 



A=A, 



d[3'(X 3 ) _ (3'jXj 
dX k 



At — A,- 



A=A, 



dX k 



X k — Xj 



dX k 
Consequently, 

(3>{X k ) KP'iXj) 8X k dX k 
Since (|8.4p implies that 

dp n (Xj) _ d 2 i3 n+l _ dp n (x k ) 



dfUXj 
dX k 



+ (5' n {X 3 



df3 n (Xj 



a=a 0'(X k )(Xj - X k ) 8X k 



dX k 



dXjdXi 



dXi 



we eventually find that (|9.10p holds for the case of j 7^ k. 

The case of j = k can be treated in much the same way, and turns out to be simpler. The 
left hand side of (|9.10p can be calculated as 

d !3' n (X k ) l^(X k ) lf3' n (X k )(3"(X k ) 



dX k f3'(X k ) 2 P'(X k ) 2 /3'(A fc ) 2 ' 
The derivatives in the parentheses on the right hand side of (|9.10p can be expressed as 



df3\X) 
dX k 



A=A fc 



df3'{X k 
dX k 



0'{X k 



:P"(X k ) 
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dX k 



A=A fc UA k 



— 7^ Pn(Afc) - -^Pnl^J- 



Thus (|9.10p turns out to hold in this case, too. 

We have thus confirmed that (|9.9p does satisfy (|9.8p . Note that (|9.9p corresponds to the 
correction terms 



n+1 ~ ^ Wi) 



for -ff n +i by the linear relation ()9.7j) . because /3^(A)'s and i?^(A)'s are linearly related with the 
same coefficients as /3 n (A)'s and i? n (A)'s. These results can be summarized as follows. 

Theorem 3. Equations of motion (|9.2p and (|9.3p can be cast into the Hamiltonian form 
The Hamiltonians are given by 

10 Examples 

We illustrate the results of the preceding section for the cases of g = 1,2,3. For notational 
simplicity, we set ^2g+i = 0. Consequently, c-i(t) disappears from various formulas. 

10.1 g = 1 

The case of g = 1 corresponds to the first Painleve equation itself. There is no higher flow (other 
than the excluded exceptional time £3). In this case, everything can be presented explicitly as 
follows. 

1) /3(A) is linear and a(X) does not depend on A: 

/3(A) = A + a(X) = ai, 

R „ u (3i sX 1 
Pi=Ri = -, ati = — = -ju x . 

2) The Darboux coordinates Ai,/ii are given by 

u 1 

Ai = -Pi = --, m = a(\i) = --u x . 

3) h(X) is a cubic polynomial of the form 

h(X) = J (A) + Ji, I (A)=A 3 + xA. 

4) The Hamiltonian J\i is equal to ifi = ii. As a function of the Darboux coordinates, Ii 
can be expressed as 

h = fj,l - Xf - xXi, 

which coincides with the well known Hamiltonian of the first Painleve equation. 
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10.2 g = 2 

This case corresponds to the two-dimensional 'degenerate Gamier system' studied by Kimura [33] 
and Shimomura [34]. A higher flow with time variables t% now enters the game. This variable £3 
shows up in the description of relevant quantities through 

c 2 (t) = 2*3- 

For instance, the linear relations between /3 n 's and Rn now read 

Pi = Ri, fo = R 2 + c 2 {t). 
Though slightly more complicated than the previous case, this case, too, can be treated explicitly. 

1) /3(A) and a(X) are quadratic and linear, respectively: 

/3(A) = A 2 + P1X + P2, a(\) = a 1 \ + a 2 . 

2) The Darboux coordinates Ai, A2,/ii,M2 axe defined as 

/3(A) = (A - Ai)(A - A 2 ), fii = a(X 1 ), /i 2 = a(A 2 ). 

3) h(X) is a quintic polynomial of the form 

h(X) = J„(A) + JiA + J 2; / (A) = A 5 + 2c 2 (t)X + xX. 

4) We still have the simple relations H\ = I\ and H 2 = I 2 between H n 's and I n 's. They are 
redefined as functions of the Darboux coordinates by the linear equations 

hXi +I 2 = nl- 7 (Ai), / 2 A 2 + I 2 = n 2 2 - I (A 2 ). 
More explicitly, 

j _ Mi ~ Jq(Ai) /i| -/q(A 2 ) 7 _ /if - Jq(Ai) a /x| -/q(A 2 ) a 
Ai — A2 A2 — A2 Ai — A2 A2 — A2 

The correct Hamiltonians K\ , K 2 are given by 

Mi M2 



^i=/i, K 2 = I 2 



Ai — Ao Ao — Ai 



Remark 10. Kimura and Shimomura studied this system as isomonodromic deformations of 
a second order scalar ODE rather than the 2x2 matrix system. In the present setting, their 
scalar ODE corresponds to the equation 

+»<*)§* +»<A»-0 

that can be obtained from (13. 7|) by eliminating the second component of ijj. The coefficients Pi(X) 
and p 2 (A) are given by 

Pi(A) = ^y, P2(A) = -h{X)-a {X) + a— — -. 

Remark 11. Actually, Kimura and Shimomura considered two Hamiltonian forms for their de- 
generate Gamier system. One of them is defined by the aforementioned Hamiltonians K%, K 2 . 
The other one is derived therefrom by a canonical transformation, and has polynomial Hamil- 
tonians. 
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10.3 g = 3 

The case of g = 3 is more complicated than the preceding two cases. Here we have two higher 
flows with time variables £3, £5. They are connected with C2(t) and 03(f) as 

5 3 
02(f) = 2* 3 ' C3 ^) = 2 tb ' 

h(X) is a sextic polynomial of the form 

h{\) = / (A)+/ 1 A 2 + / 2 A + / 3 . 
A new feature of this case is the structure of Iq(\): 

J (A) = A 7 + 2c 2 (f)A 5 + 2c 3 (t)A 4 + {c 2 {tf + x)X 3 . 

Note that the coefficient of A 3 is now a quadratic polynomial of the time variables. Of course, if 
we consider the general case (|6.9|) . this is a rather common situation; the first two cases (g = 1 
and g = 2) are exceptional. 

11 Conclusion 

We have thus elucidated the Hamiltonian structure of the PI hierarchy for both the Lax equa- 
tions and the equations of motion in the spectral Darboux coordinates. Though the extra 
terms U' n+l (\) in the Lax equations give rise to extra terms in the equations of motion for 
the Darboux coordinates, these terms eventually boil down (somewhat miraculously) to the 
correction terms /S.H n+ \ in the Hamiltonian. 

The correction terms AH n+ i are identified by brute force calculations. It is highly desirable 
to derive this result in a more systematic way. As regards the Gamier system, such a systematic 
explanation is implicit in the work of the Montreal group |28[ [29], and presented (in a more 
general form) by Dubrovin and Mazzocco [31]. Let us recall its essence. 

As mentioned in Introduction, the Gamier system is equivalent to the 2x2 Schlesinger system. 
The L-matrix of the Schlesinger system is a 2 x 2 matrix of rational functions of the form 

N . 

^(A) = £ — 



1 A - h 



The matrix Aj takes values in a two-dimensional coadjoint orbit of sl(2,C). This orbit, as 
a symplectic leaf, carries special Darboux coordinates £j, rjj. Aj can be thereby written as 



2 



1 ( dm q 



Aj A '/;' • "jV -to* 

where Oj is a constant that determines the orbit, and may be interpreted as a monodromy 
exponent at the regular singular point at A = tj. The Lax equations can be converted to a 
Hamiltonian system in these Darboux coordinates rjj. Since the spectral Darboux coor- 
dinates Aj, fij are connected with these Darboux coordinates by a time- dependent canonical 
transformation, the Hamiltonians in the latter coordinates have extra terms. 

Unfortunately, this beautiful explanation of extra terms does not literally apply to the present 
setting. The relevant Lie algebra for this case is not sl(2, C) but its loop algebra sl(2, C)[A, A" 1 ], 
coadjoint orbits of which are more complicated. 

A similar idea, however, can be found in the recent work of Mazzocco and Mo |41j on an 
isomonodromic hierarchy related to the second Painleve equation. They start from the Lie- 
Poisson structure of a loop algebra, and convert the Hamiltonian structure on a coadjoint orbit 
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to a Hamiltonian system in Darboux coordinates by a time-dependent canonical transformation. 
It will be interesting to reconsider the present setting from that point of view. 

Another remarkable aspect of the work of Mazzocco and Mo is that they present another 
set of Darboux coordinates alongside the spectral Darboux coordinates. Unlike the spectral 
Darboux coordinates, these coordinates are rational functions of the dynamical variables in the 
Lax equations; this is a desirable property in view of the Painleve property of the system. 

Actually, borrowing their idea, we can find a similar set of Darboux coordinates Q n , V n for 
the PI hierarchy as 

where p^ stands for the fc-th power sum 

3=1 

(Note that p^s and (3^s are related by the generating functional relation 



Pk , /3(A) 
log- 



oo 

2 kX k ~ ~ 6 \9 ' 

k=l 

so that pk's may be thought of as polynomial functions of /%'s.) One can prove, in the same 
way as the case of Mazzocco and Mo, that Q n , V n do satisfy the canonical Poisson relations. In 
fact, they turn out to satisfy the stronger relation 

9 9 

d\j A dfij = dQn A dP n , 

j=l n=l 

which implies that Q n , V n are connected with the spectral Darboux coordinates Xj, pj by 
a time-independent canonical transformation. In particular, the canonical transformation yields 
no correction term to the transformed Hamiltonians. Namely, K n+ is persist to be correct 
Hamiltonians in the new coordinates Q n , V n as well. Moreover, employing the Lagrange in- 
terpolation formula, one can see that K n+ i's are polynomials in these coordinates (and the 
time variables). We thus obtain a generalization of Kimura's polynomial Hamiltonians for the 
degenerate Gamier systems [33]. 
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